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Abstract 

We analyse the properties of generic models based on an SU (3) family symmetry pro- 
viding a full description of quark charged lepton and neutrino masses and mixing angles. 
We show that a precise fit of the resulting fermion textures is consistent with CP being 
spontaneously broken in the flavour sector. The CP violating phases are determined by the 
scalar potential and we discuss how symmetries readily lead to a maximal phase controlling 
CP violation in the quark sector. In a specific model the CP violation to be expected in 
the neutrino sector is related to that in the quark sector and we determine this relation 
for two viable models. In addition to giving rise to the observed structure of quark and 
lepton masses this class of model solves both the CP and flavour problems normally as- 
sociated with supersymmetric models. The flavour structure of the soft supersymmetry 
breaking terms is controlled by the family symmetry and we analyse some of the related 
phenomenological implications. 



1 Introduction 

The gauge sector in the Standard Model (SM) is completely fixed by gauge invariance once 
we specify the particle content and its quantum numbers. This sector has been thoroughly 
tested with great precision at the collider experiments and the SM has been confirmed as the 
correct gauge theory of strong and electroweak interactions up to energies of order 100 GeV. 
In contrast, the flavour sector of the SM remains the big unknown in high energy physics. In 
the SM flavour (Yukawa) couplings are not determined by symmetry and so they are merely 
unrelated parameters to be determined by experiment. In our search for a more fundamental 
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theory we need to improve our understanding of the flavour physics. In analogy with the gauge 
sector, one might expect a spontaneously broken family symmetry determines the different 
flavour parameters. This approach has already been tried several times with reasonable degrees 
of success pQ although so far it is difficult to choose one of the various options on offer. The 
appearance of CP violation in the SM is equally mysterious. We do not know why there are 
complex couplings in the SM giving rise to a violation of the CP symmetry in neutral meson 
systems. In fact, these two problems are deeply related in the SM as the Yukawa couplings are 
the only source for both the flavour structures and the CP violation phenomena, although this 
need not be so in extensions of the SM where new sources of CP violation different from the 
Yukawa couplings are present. 

In this paper, we address both problems together and build a supersymmetric flavour theory 
which determines all the flavour structure of the theory as well as the different phases breaking 
the CP symmetry. This is particularly relevant in the supersymmetric context where the lack 
of understanding of flavour and CP is especially severe. The minimal supersymmetric extension 
of the SM (MSSM) has many new scalar states which introduce the possibility of new sources of 
flavour changing neutral currents (FCNC) beyond the usual Yukawa matrices giving rise the so- 
called "SUSY flavour problem". Similarly there are new CP violating phases that, unless small, 
give unacceptably large contributions to CP violating observables - the "SUSY CP problem" 
The simplest solution to these problems consists of arbitrarily assuming that all the new 
flavour structures are proportional to the identity and that all new phases vanish. However, this 
approach cannot be justified without some underlying symmetry reason (so far unspecified), 
made more difficult to understand given that the Yukawa couplings of the SM do not share 
these features. 

We will study both CP violation and the origin of flavour structures in the framework of a 
supersymmetric SU(3) family symmetry model. This kind of model has been shown [3J HJ |3j to 
provide a correct texture for the Yukawa matrices in agreement with an accurate phenomeno- 
logical fit 0. However, the improvement on the determination of CKM mixing angles and CP 
asymmetries requires the presence of phases in the elements of the Yukawa matrices. In our 
analysis we assume that CP is an exact symmetry 4 of the theory of unbroken flavour SU(3) 
and is only spontaneously broken by complex vacuum expectation values of the flavon fields that 
determine the Yukawa structures jH]. As we will discuss this approach has the great advantage 
of naturally explaining the SUSY CP problem. Moreover the CP violating phases are deter- 
mined at the stage of flavour symmetry breaking and we will discuss how maximal CP violation 
naturally results. 

In Section 12 we present the general features shared by SU(3) flavour models which reproduce 
the observed quark textures. In Section El we show that spontaneous CP violation in the flavour 
sector is able to describe quark masses and mixing angles while generating the observed CP 
violation. In Section 0] we show how the same structure can reproduce the correct masses and 
mixing angles in the leptonic sector. We give two different examples of SU(3) models which 
generate acceptable neutrino masses and mixings and determine the prediction for the different 
CP violation observables in the lepton sector. In Section El we show that these models solve 
both the supersymmetric flavour problem and the supersymmetric CP problem and we point 
out other observables where signals of a supersymmetric SU(3) flavour model may show up in 

4 In string theory CP invariance can be a symmetry of the 4D compactified theory corresponding to an 
unbroken discrete subgroup of the higher dimensional Lorentz group |7j 
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future experiments. Then, in section El we analyse the vacuum alignment giving rise to the quark 
and lepton flavour structures and finally we present our conclusions. 



2 SU(3) flavour symmetries and Yukawa textures 

One reason why it is difficult to construct a theory of flavour is that measurement in the quark 
sector of the eigenvalues (quark masses) and the CKM mixing matrix is all the information we 
can extract about the full quark Yukawa matrices using SM interactions alone. Unfortunately, 
this is not enough to determine the full structure of these matrices. Therefore there remains a 
lot of freedom corresponding to widely different Yukawa textures, ranging from democratic to 
strongly hierarchical matrices which, in the SM, give rise to identical phenomenology. However, 
in most extensions of the Standard Model, and in particular in super symmetric extensions, the 
new non SM interactions can distinguish between the different options. 

Before the discovery of new supersymmetric interactions or any other new interactions able 
further to probe the flavour sector we must rely on simplifying assumptions to try to disentangle 
the complicated structure of masses and mixing angles. A particularly interesting approach is 
to look for texture zeros [9j|, small entries in the Yukawa textures which lead to viable relations 
among the eigenvalues and mixing angles. An example is the well known Gatto-Sartori-Tonin 
relation JU] and its phenomenological success may indicate some underlying dynamics capable 
of generating those zeros. Given the success of the GST relation we consider it reasonable to look 
for textures with this property. Moreover the smallness of left handed mixing angles suggests 
that the non-diagonal elements are smaller than the diagonal elements. A recent phenomeno- 
logical analysis |E] under these assumptions shows that the following symmetric textures give 
an excellent fit to quark masses and mixing angles, 
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with e ~ 0.15, e ~ 0.05, b = 1.5, a = 1.3, c = 0.4 and a', b', d are poorly fixed from experimental 
data. Note that with the assumption of small off-diagonal elements the data determines those 



entries above the diagonal (in our convention Qid c Y^ 



In this paper we adopt this basic structure as our starting point. As shown in [3J |3] these 
structures can be successfully reproduced from a spontaneously broken SU(3) family symmetry. 
Under this symmetry all left handed fermions (ipi and ipf) are triplets. To allow for the spon- 
taneous symmetry breaking of SU(3) it is necessary to add several new scalar fields which are 
either triplets (6*3, 6*23, 9 2 ) or antitriplets (9 3 , #23)- We assume that SU(3)fi is broken in two 
steps. The first step occurs when 9 3 gets a large vev breaking SU(3) to SU(2). Subsequently a 
smaller vev of #23 breaks the remaining symmetry. After this breaking we obtain the effective 
Yukawa couplings at low energies through the Froggatt-Nielsen mechanism ^T] integrating out 
heavy fields. In fact, the large third generation Yukawa couplings require a 83 (and 83) vev of 
the order of the mediator scale, Mf, while 823/Mf (and 823/Mf) have vevs of order e in the up 
sector and e in the down sector with different mediator scales in both sectors. Through vacuum 
alignment the D-terms ensure that the fields #23 and #23 get equal vevs in the second and third 
components. This is necessary to generate the Yukawa structure of Eq. (JTJ). Following 0, we 
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assume that 8 3 and 6*3 transform as 3 © 1 under SU(2) R necessary to differentiate the up and 
down Yukawas. The SU(2) R arises naturally in the context of an underlying SO (10) grand 
unified theory and with this it is possible simultaneously to describe quark and lepton Yukawa 
matrices. The main ingredient to reach this goal, is a new field (or composite of fields) S, which 
is a 45 of 577(5) with a vev along the (B — L + k,T^) direction with k = 2, corresponding to 
weak hypercharge. This generates the usual Georgi-Jarlskog factors that correct the difference 
between down quark and charged lepton Yukawas. 



2.1 The superpotential 

The SU(3) and 5*0 (fO) symmetries are not enough to reproduce the textures in Eq. and 
we must impose some additional symmetries to forbid unwanted terms in the effective superpo- 
tential. The choice of these symmetries is not uniquely defined (two examples are presented 
in Section 4). Nevertheless the quark Yukawa textures consistent with the requirement of 
SO (10) x SU(3) unification strongly constrain the superpotential describing the Dirac Yukawa 
couplings. In an SO(10) x SU(3) model, third generation Yukawa couplings are generated by 
8181 while couplings in the 2-3 block of the Yukawa matrix are always given by ^23^23^- We 
have two options for the couplings in the first row and column of the Yukawa matrix. They can 
be given by e lkl 823^83,1823 (82383) or e tkl 8 23^83, 18 23 (#3 #23) , although the requirement of correct 
Majorana textures favours the first option in all the cases analysed, as shown in Section 0] With 
this the basic structure of the Yukawa superpotential is given by 5 

w Y = [eie{ + ei 3 ei 3 E + e ik % 3)k 8 3)l ei 3 (# 23 ^) + 

e ljk 8 2 3,k {82383') 2 + e ij %, k (823T3) (8238V3) + • • 

This structure is quite general for the different 577(3) models we can build, although in particular 
examples we may have different subdominant contributions that can play a relevant role in some 
observables. As the structure of Eq. (j2J) fixes the main features of mixings, eigenvalues and CP 
phases we can make some general statements about the phenomenology of such models. 

On minimisation of the scalar potential the SU(3) symmetry aligns the vevs giving the form 

\ a _ ( n\ ( a 3 ^ 
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As we will see in section |H1 these vevs are of a form such that, 

— = — = e 1/4 — = e — = e (4) 

where, to fit the masses and mixing angles, e ~ 0.15 and e ~ 0.05 at the symmetry breaking scale 
that we take approximately equal to Mq UT . Notice that, although a 3 /M = e 1//4 is quite large as 

5 The coefficients of O(l) of the various terms have been suppressed for clarity 
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an expansion parameter, due to the SU(3) symmetry it only enters in the low energy parameters 
as (9 3 /M) 2 = e 1 ! 2 or {6 23 6 3 ) / M 2 (see for instance Eq. (J2J)). Then, the most important parameter 
for the hierarchy of fermion masses is #23 which is quite small and the series expansion converges 
rapidly. So, with these vevs we have the following Yukawa textures, 



( 



V ... 



s 3 e i6 X 1 e 3 e^ f+ ^X 2 \ 



\a 



E 

7J2 



£ 2 e iPs 



\a 



E 

TJ2 



J 



\a{\ 2 

M y 



(5) 



where X\ and X2 account for the different contributions to the (1,2) and (1,3) elements from 
the last three terms in Eq. (J2J) and 5 u,d = (2a Ut d + P3 + /3' 2 ). As mentioned above the vev of E 
changes for different fermion species. The vev a 3 and the mediator masses, Mf, are different for 
fermions in the up or down sector while for the up quarks and neutrinos \ = 0. 

In the same way, the supergravity Kahler potential receives new contributions after SU(3) fi 
breaking. As explained in the introduction, we assume that CP is an exact symmetry of unbro- 
ken SU(3) fi and it is only spontaneously broken by the complex flavon vevs. Thus, at M P i anck 
both SU(3)fi and CP are exact symmetries of the theory and both the superpotential and the 
Kahler potential are invariant under these symmetries. Furthermore, we use a Giudice-Masiero 
mechanism to generate a /x-term of the correct magnitude so we have in the Kahler potential a 
real coupling mixing the two Higgs doublets. In this class of theory the hidden sector F-term 
breaking SUSY is also real. Therefore, after SUSY breaking, we obtain a real /i-term. In Section 
|5]we show that the corrections after SU(3)fi breaking do not modify this conclusion. 



2.2 The Kahler potential 

After the flavour symmetry is spontaneously broken we obtain the Yukawa structure of Eq. (jSJ). 
In the same way we can calculate the effective Kahler potential that will be a general nonrenor- 
malizable real function invariant under all the symmetries of the theory coupling the fields ip\ipj 
and flavon fields. In first place we must notice that a term is clearly invariant under gauge, 
flavour and global symmetries and hence give rise to a family universal contribution. However, 
SU(3) breaking terms give rise to important corrections ^21 IT3| . Any combination of flavon 
fields that we add must be also invariant under all the symmetries, although the SU(3) indices 
can also be contracted with quark and lepton fields. In fact, it is interesting to notice that, 
due to the possibility of introducing the hermitian conjugate fields in the Kahler potential, new 
terms are allowed with different suppression factors from the terms that appear in the Yukawa 
couplings of the superpotential. This has important effects both in the fermion mass matrix 
and in the scalar mass matrix. To include the latter we need to know the mechanism of super- 
symmetry breaking. We assume that it is broken in the hidden sector through a field X with 
non-vanishing vacuum expectation value for its F-term. Including this field the leading terms 
in the Kahler potential for the matter fields are given by 

K = 4^ ^(co + d XXt) + ^[0*4(01 + diXYt) + 44(c 2 + d 2 XX^)] + 

+^^ k % tk 9 23 , i y^ mn \ m 9 23>n )(c 3 + d 3 XX^)j (6) 
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where ip represents the SU{2) doublet or the up and down singlets and accordingly Mf = 
Ml, M u , Adj. Then Cj,cfj are 0(1) coefficients in the different terms. After SU(3)fi breaking 
we obtain wave function and soft mass corrections because the fields acquire non-zero vevs and 
F-terms. 



2.2.1 Wave function effects 

Due to the terms independent of X, the matter fields do not have canonical wave functions in 
the basis 1/1%, ipj [Hj. Rather, in the case of the SU(2) doublets, they are given by 

c + c 3 a 2 b 2 23 /M L 

(4 4 4)\ c Q + c 2 bl 3 /Ml c 2 bl 3 e-^/Mt 

c 2 b 2 23 e^/M 2 c + Cl (af + a d 3 2 )/Ml 

= (4 4 4)^n\ ^ 2 I (?) 
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where the messenger mass Ml corresponds to heavy SU(2)l doublet fields and the coefficients 
Ci are of 0(1). There is an analogous expression for the wave functions of the conjugate fields ip c 
with M L replaced by M Ujd the mass of the u, d SU(2) singlet messenger masses and a 3 2 + a 3 2 
replaced by 03 2 or a 3 2 for the up and down sectors respectively. As discussed in Ml £3 M u > 
Md is necessary if the expansion parameter for the up quarks is to be smaller than that for the 
down quarks and the Majorana mass matrix is to have a phenomenologically acceptable form. 

The effect of these wave function corrections on the Yukawa couplings may be determined 
by redefining the matter fields to obtain canonical kinetic terms, ip' = Nip. Consider first the 
contribution to the observed mixing angles from the ip wave function effects. From Eq. (j3J) we 

have t| = e 3 . In this case we have ip' = Nip = PJX^Pvp where 
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This is common to the up and the down sectors. As a result, if X were the unit matrix, the 
effect of R would be to give a common contribution, R\ to the matrices diagonalising the up 
and down quark mass matrices, V u and Vd- This can be absorbed in a redefinition of the LH 
multiplets and clearly does not generate a CKM element. Thus the leading contribution from 
these terms to the CKM matrix comes from the departure of X from the unit matrix and is thus 
0(e) suppressed and occurs at 0(e 4 ). Thus the mixing introduced in the (2,3) sector through 
the Kahler terms involving the fields ip is much smaller than that coming from the superpotential 
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|15j . Turning to the Kahler potential of the fields ip c sector, due to the ordering M u > the 
mixing introduced is different for the up and the down sectors and is larger in the down sector. 
Due to this there is a mixing of 0(e 2 ) introduced in the ip c sector from the Kahler potential 
|15j . This mixing in the right handed sector can also modify the CKM mixing, however this 
effect is strongly suppressed. If the matrix X relevant for the ifi c sector is the unit matrix the 
effect vanishes. Including the effect of the higher order terms responsible for making X different 
from the unit matrix we find a suppression by the factor (m c /m.&)(— r 2 + r 2 ) times the mixing 
diagonalising the right handed sector 0(e). Thus, this results in a contribution again of 0(e A ) 
and much smaller than that coming from the superpotential [To] . 

One may readily check that the effects coming from the Kahler potential also do not affect 
the MNS matrix to leading order. 



2.2.2 Soft mass corrections 



The Kahler potential Eq. © also gives corrections to the soft SUSY breaking scalar masses gener- 
ated by the vev for the F-term of the X field. Redefining the matter fields to obtain the canonical 
kinetic terms we obtain the soft breaking mass matrices as, M?- = F a F b (d log Kij) / (d(p a d(f) b ). 
In general the coefficients q will not be proportional to di and therefore diagonalising and 
normalising the kinetic terms will not diagonalise the sfermion mass matrices. Then in this 
model, suppressing factors of order 1 due to the non-proportionality of q and di, we have the 
contribution to the scalar masses of the form 
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where m = \F X \ 2 / M P i anck and Mf = M L , M u , M d for the doublet squarks and the u, d singlet 
squarks respectively. 



3 Phenomenological analysis of phases and CP violation 
in the quark sector 

The general structure of the Yukawa matrix and the soft breaking mass matrices in the SU (3) 
family symmetry models is given by Eqs. (0) and (jllj) . The parameters are fixed mainly by the 
quark masses and mixing angles and in the fit to data of the quark textures the presence of 
phases plays an important role. Furthermore CP violation has been measured experimentally 
in both the neutral kaon and neutral B sectors and we must reproduce these measurements. 

To address this we turn to a detailed discussion of the phenomenological fit. It is convenient 
to use the Wolfenstein parameterisation for the CKM matrix with parameters A, p and rj. We 
have 

§^ = A V / 7T?, |^ = -v / ( C A-p) 2 + f, \V US \ = X, (12) 
\V cb \ c x \V ts \ c x 

for p = cxp and f] = c\T], cx = 1 — A 2 /2. The experimental constraints are the semileptonic 
decays of B mesons, whose charmless channel determines the ratio |T4&/V^,|, the CP violation 
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in the K system, from which we obtain the asymmetry parameter |e#|, the measurements of 
B® s — B® s oscillations, from which the mass differences Am^ and Am# s can be obtained, and 
the CP asymmetries in various B decays, from which sin (2/3) is obtained with (3 one of the angles 
of the unitarity triangle. We have updated the 2001 fit with the latest experimental information 
6 . In this fit we have included the constraints IV^/V^I, ArriB s , Am^^ and \ck\- In the limit 
where we neglect all supersymmetric contributions to these observables, the fitted values for p 
and fj are 
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(13) 



We prefer to leave aside the constraints coming from the decay modes of B, given the present 
uncertainty regarding the decay mode B — > <&Ks- 

Let us comment on how the phases and the mixing angle sf 3 are fixed. For matrices with 
the hierarchy of the form of Eq. (jSJ), the magnitudes of the ratios V t d/V ts and V u t,/V c b depend 
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sf '■ and sVj the rotations diagonalising the down and up quark Yukawa textures respectively and 



j{j = Arg(if^) = ^g(^2k m ik m jk)- The ra ti° is the most sensitive to the presence of a 
phase in Eq. (fT4*j) due to the contribution coming from sf 3 /s^ 3 . At first approximation, the phase 
$i can be fixed by using the GST relation 
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(16) 



where we have used that s" 2 = a/^t an d sf 2 = y^S which applies for the structure of the mass 

matrix in Eq. (j3J). The ratios of masses m u /m c and md/m s can be written in terms of ratios 
^uj^d and the parameter Q, which are better determined from chiral perturbation theory JE] 
(see table 3 in jH]) and the ratio m s /m c can be obtained from the MS values at a common scale 
|17| ITS! ITT?] . This information gives us 



— = 0.228 ±0.005, 



TYl 

— = 0.054 ±0.011. 



(17) 



Comparing Eq. (|16j) with the experimental value of V us = 0.2240 ± 0.0036, we obtain $! 
± (1.3 ± 0.6) which corresponds in degrees to $i = ± (74° ± 34°). 



3 The input data used in this fit are as given in in Appendix A, Tableland in Table 1 of 
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From here, we see that non zero phases are needed to fit quark masses and CKM mixings, 
even before considering CP violation observables. Taking into account that $1 is the main phase 
contributing to the Jarlskog invariant, as we will see below, this implies that its value is in the 
upper range of this interval. 

However in order to fix appropriately the phases and the mixing angle sf 3 , we need to give 
a common solution to the system of the two equations appearing in Eq. (JT4"jl and to Eq. (fTHj) 
for the variables $2; $1 and sf 3 , and we need to use a constraint on $ 3 in order to write the 
equations just in terms of the above variables. For matrices of the form Eq. (JTJ), $! $ 2 and $3 
are given by, 

$1 « (7i u 2 - 7i 2 ) - (722 - 722) 

$2 « (7i d 3 - 723 - 7i d 2 + 72 2 )-^sin( 72 rf 3- 7 ^ 3 )-sin(^)4^ 

S 23 

$3 « (7i d 3 " 7 2 d 3 - 7r 2 + 722) 4^ 3 - 7 2 3)-^)44 (18) 

s 23 

with </7 = 7(3 — 7(2 — 723, the coefficients c{ 2 are of order 1, and as defined above 7^ = 
Arg(^ fc m{ k rrij^). Thus we can write 

$ 3 «($2-$l). (19) 

Hence the ratios Eqs. (|14|lfij) of the CKM matrix, can be expressed just in terms of two phases, 
$1 and $ 2 



V ts \ \l m a V m 



|r "' 1 - / ^ + ^-2^/^cos($2-$ 1 ; 



11 "''I - /^ + C2_ 2 C A &cos($ 2 ) (20) 



I I \/ m r . V m 

for C = c£° (with c the 0(1) coefficient in Eq. (0)), for e° = 0.15 ± 0.01 7 . Solving for the 
three equations, Eqs. (|16I20|) we have the following solution, requiring c to be real and to satisfy 
< c < 1, 

$1 =1.57 ±0.02 — ► c = 0.28 ±0.09, $ 2 = 0.8 ± 0.3 (21) 

Notice that the inclusion of the constraints in Eq. (|2U|) modifies substantially the mean and 
uncertainty of $1 from the result we obtained just from Eq. (|T5j). We have included all the 
uncertainties of the parameters involved in Eqs. (|14|) and (|16j) 8 . 

3.1 Comparison with the 577(3) texture. 

These values must be compared with the predicted values in our model from Eq. (J18)) and the 
leading order texture in Eq. (jSJ), 

$1 = 2{a u - a d ) (22) 

7 The quantities with superscript refer to the quantities appearing in the parameterisation of Eq. (25) in . 
8 We can compare these solutions with those of the fit of [2] where p, r\ and A are different: c = 0.5, $2 = —24°, 
c = 0.95, <J> 2 = -(58 ± 5°) for $1 = +90° and $1 = -90° respectively. 
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* 2 «--fsiii(7&-7£,) (23) 

s 23 

where we used the fact that, at leading order, 723 = /3 3 — 2%, jf 2 = ^ 7i3 = & d + /3 3 — 2%, 
723 = 7i2 = ^" an d 7i3 = & u + ^3- So, we obtain <E> 2 ~ 0.11 assuming that the phase 
723 — 723 = 8 — 5 U — 2x is 0(1). Therefore, from the leading order texture in Eq. (J3J we can 
only be compatible at the 3a level with the required value from the SM fit which corresponds 
to $ 2 = 0.771™. 

However, there are several other contributions which can affect this result. In particular 
there are subleading contributions in the e expansion to the different elements of the Yukawa 
texture and also there may be significant supersymmetric contributions to SM processes which 
can affect the fit to $ 2 . 



3.2 Subleading contributions 

In the case of Model II analysed later j^J, we have an additional contribution to Y 2 i such that, 

K? 3 = ^32 = ? e«^ + c W (24) 

l a 3l 

with c an unknown coefficient order 1. Then we have, 

$2 - c ^ ^ 5 ( 1 - -^-) sin(/3 3 - B S ) + O(^) (25) 

Therefore, in this case we have a contribution of order yl ~ 0.4 which can accommodate the 
result of the fit within la of $ 2 - 



3.3 SUSY contributions 

Another possibility is that there are significant SUSY contributions to SM processes that have 
not been taken into account in the SM fit and can play an important role in determining the 
value of $ 2 . A full computation of the complete SUSY contributions is beyond the scope of this 
paper, but we can estimate the size of the new SUSY contributions and fit a modified standard 
model triangle. 

Using Eqs. (0) and (jllj) we can obtain the sfermion mass matrices in the so called SCKM 
basis with the correct phase convention in the CKM mixing matrix. The SCKM basis is the basis 
where we rotate the full supermultiplet, both quarks and squarks, simultaneously to the basis 
where quark Yukawas are diagonal and real. Then, we must also ensure that the phases in the 
CKM matrix are located in the (1, 3) and (3, 1) elements up to order e 4 . This is done explicitly 
in Section El and here we simply use these results. To analyze flavour violating constraints at the 
electroweak scale, the model independent Mass-Insertion (MI) approximation is advantageous 
[2U] . A Mass Insertion is defined as 8^ = A^/m^; where A{- is the flavour- violating off-diagonal 

entry appearing in the / = (u, d, I) sfermion mass matrices in the SCKM basis and m?~ is the 
average sfermion mass. In addition, the MI are further sub-divided into LL/LR/RL/RR types, 
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labeled by the chirality of the corresponding SM fermions. The different Mass Insertions at 
the electroweak scale relevant for ek are, 



|Im(^)f 2 | = 6.8 x l(r 4 sin$i (26) 



|Im(4 fi ) 12 (^ L ) 12 |=2xlO- 4 sin$ ] 



These values are only order of magnitude estimates representing the maximum possible value 
for these entries and they can be reduced by factors order one in the redefinition of the fields 
to get the canonical kinetic terms as explained in Section This must be compared with the 
phenomenological MI bounds [20] . 



\M< R )i2\ < 3-2 x 10 



-3 

|Im(^ R ) 12 (^ L ) 12 |<2.2xlO- 4 (27) 

These bounds correspond to average squark masses of 500 GeV and they scale as rriq (GeV)/500 
for different squark masses. Then, for rrig ~ 500 GeV, it is clear that we can assume a sizeable 
SUSY contribution to Ek that can reach a 20% of the experimental value. For larger squark 
masses a 20% contributions is still possible if we consider simultaneously (^r)i2 and (5^)i2- 

To take into account of a possible SUSY contribution up to a 20% of the experimental value, 
we enlarge the error on ek to 20% of the mean value. The resulting fit is given in Figure ^ giving 

p = 0.212±g52, ?7 = 0.348^2, (28) 

which in turns imply different values for the parameters c, $ 2 : 

$i = 1.58 ± 0.02, c = 0.297±g;g||, $ 2 = 0.830±^ (29) 

Although here the errors on <£> 2 have increased, so the mean, which again would be compatible 
with the result of Eq. ()23j) only at 3a level. This shows that the favoured solution prefers a 
SUSY contribution of opposite sign to the SM contribution. 

However it is instructive to determine the implication of a supersymmetric contributions to 
€k which has positive sign with respect to the SM contribution, such that the SM value of €k is 
only 80% the mean of the experimental value of ek- Fixing this a fit to the revised data gives 

n - 228+ ' 057 ft - S06+ - 050 

r — U ' ZZO -0.059> 7 — U - OUO -0.052 

$i = 1.57 ± 0.02, c = 0.3361^1, $ 2 = 0.599^J£ (30) 

Although in this case the quality of the fit is somewhat smaller than in the previous fit, now the 
model I is compatible within 2cr of the result of the fit for $ 2 -Eq. (|30p. while the model II is 
compatible at la of the fitted value. 

3.4 The Jarlskog invariant 

It is instructive also to calculate the Jarlskog invariant in the CKM matrix which can be com- 
puted from 

Jcp = lm{V 3 * 2 V 2 \V 22 V 31 } (31) 
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where V is the CKM matrix, which depends on the sines of the mixing angles, sin^(0y) = s^-. 
In this case we have 



(32) 



r// iff -.// 

I f mi, mi, i if, rni 9 rni, 9 i i f mJ.m , 

H — f moo H — 32 . 22 \min-\ — 3 f 

J „ m 22 J „ m 33 J _ m 33 

\ m 22\ \ m 33\ \ m 33\ 

which to leading order are given by 

s{ 2 = 0(e), 4s = 0(e 2 ), 4 = 0(e 3 ), (33) 
with e = e for / = c/-quarks and e = e for / = w-quarks and e ^> e. Thus Jcp acquires the form 



J CP w Im{s 



23 



Q d Q U — Q d Q d P^ 2 — e" Q d 



*2 3 e-^ = s^e- l ^3-7 23 ; _ s « 3 (34) 

with real. 

In a particular parameterisation of the CKM mixing matrix we have a single unremovable 
CP violating phase which is responsible for CP violation through this invariant. In the PDG 
"standard" parameterisation [21] the Jarlskog invariant is Jcp = S12C12S23C23S13C13 sin 5, with 
Cij,Sij the cosine or sine of a rotation in the plane. Using the measured values of the 

mixing angles at the electroweak scale we obtain Jqp = 3.3 x 1CT 5 sin<5 ~ (2.8 ± 0.4) x 10 -5 , 
with 5 = 1.02 ± 0.22. In our model the Jarlskog invariant in Eq. (|3*4*|l must reproduce this value 
up to a correction in the phase due to the relatively small SUSY contributions to the unitarity 
triangle fit (see Figure^). We must remember that $1 is nearly maximal, while $2 is at most 
order a/1 or smaller. Using this together with the fitted values of the coefficients from Eq. (0) 
this implies that the Jarlskog invariant is dominated by the first contribution in Eq. (|34j) . If we 
use the values of the angles at Mqut we obtain, 

j cp\ Mgut ~ a 2 b b'e 5 esm^ 1 ~ 1.44 x 10" 5 sin$i (35) 

where we used that b' = b as the (1,2) entries are not affected by the vev of S breaking SU (2)r. 
However, this is the value of the invariant at Mqut and we need to evolve it to Mw to compare 
with the experimental results. We can use the RGE equations for the Jarlskog invariant in 
and we obtain, 

^cp| Mw - X" 2 ^cp| MG[/T ^2.9xlO- 5 sin$ 1 (36) 

where x = {M GUT / M Z )~ Y ? ~ 0.7 [23 Ej- By comparison with the experimental result for 
the Jarlskog invariant we again conclude that $1 must be near maximal. 



4 Neutrino masses, mixings and CP violation 

As we have seen in the previous Sections the quark Yukawa couplings fix many of the features of 
our SU(3) fi model, but we still have some freedom to build the right- handed Majorana neutrino 
mass matrices. To do this we must take into account that this matrix violates lepton number by 
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Table 1: Transformation of the matter superfields under the family symmetries. 



two units. As discussed in Appendix C we assume this breaking comes from a sneutrino triplet, 
A, and antitriplet, A with vevs A = (0, 0, c) and A = (0, 0, c). Then the dominant Majorana mass 
comes from 

W M = i#Aty?A' + . . . (37) 

Any other term in the superpotential will have to include this set of fields because A is the only 
term violating lepton number by 1 unit. Therefore all other terms are obtained by combining 
this term with a neutral combination of fields and reordering the SU(3) indices. The allowed 
terms will depend on the global symmetries that must reproduce the correct quark and charged 
lepton Yukawa textures and at the same time give rise to a suitable Majorana mass matrix. 
As mentioned in Section I the choice here is not unique and we have several possibilities. Here 
we will give two concrete examples which reproduce the correct Yukawa textures and give an 
acceptable neutrino mass matrix. Since CP is spontaneously broken only in the flavour sector 
all the phases in the lepton sector are a function of the flavon phases. Consequently there are 
predictions for the CP violating phases in neutrino physics. We will illustrate this in the two 
examples given below. 



4.1 Model I 

In this example we need two additional global symmetries to obtain the correct Yukawa textures. 
These include a continuous R-symmetry which can play the role of the Peccei-Quinn symmetry 
24 to solve the strong CP problem and a discrete symmetry Z 3 <g> Z 5 , under which the fields 
transform as shown in Table ^ 

4.1.1 Dirac structure 

Using these charges we can build the effective Yukawa couplings as, 

W Y = [^3 + + (^23,^3,^23 + ^23,* V23) (#23#3~) 

e ij % 3>k {02 3 9~ 3 ) 2 + e ij %, k (9 23 T 3 ) (9 23 9^) + (^*] (38) 

and this reproduces the structure in Eq. (JSJ), apart from the last term that will give a subdomi- 
nant contribution to the (2, 3) and (3, 2) elements of the Yukawa matrix. The neutrino Yukawa 
matrix is given by Eq. (|3*5jl with Y, v = 0, 



e 3 e i6u X 1 e 3 e^ u+ ^X 2 . „,, 

Y V = I ... v£ 2-75 e iia u l|3l ( 39 ) 



1 



Ml 



13 



with 7] an order 1 coefficient. 



4.1.2 Majorana structure 

In a similar way we can calculate the structure of the Majorana matrix. As explained above 
we simply have to look for the neutral combinations under the extra symmetries. We have two 
neutral combinations of fields that will be relevant: (#23^3) (9 3 9 3 ) and (#23^3 ) ($3^3) 9zi@3j@23k 
(the combination (9 3 9 3 ) 9 3i 9 3 j9 23k om y contributes to the element (1,3) and has a small impact 
on neutrino mixings). So, we obtain the different entries in the Majorana matrix from, 



W: - - ^ J 



c„i.c 

x'x j + eiei 3 (a#3~) 2 (9 3 r 3 ) 2 + Q\A (a^) 2 («)° (40) 



M 



M 

+ e lk % 3 J 3tl 9i (\9~ 3 ) 2 (9 3 T 3 ) 6 + e ik % 3 ^ 3 ^ 23 {\W 3 ) 2 («) {9 23 T 3 )' 



and assuming b 23 /M u = £ and a 3 /M u = e , this gives a Majorana matrix of the form, 



Xe 5 - 75 e iS ' e 5 e i(S'+3a u -2/3 s ) 



c 2 



Mr = ( ... ^ 2£ 5.5 e i(8a u ) ^2.75^(4^) j _ ( 41 ) 

with 5' = (6a u + 2f3 3 + (3' 2 ) and A, 771, r] 2 , Oil) coefficients. 

Using Eqs ()39|) and ()41|) . the form of the low energy neutrino mass matrix is determined in 
Appendix B, c.f. Eq (jHHj) . together with the corresponding neutrino masses and mixing angles. 

For a reasonable choice of the parameters we can easily obtain the observed masses and 
mixing angles. Choosing X\ = 1.4, X 2 = 2, rj = 1.5, A = 0.34, r] 2 = 0.75 and 771 = 1/3 and 
neglecting charged lepton mixings, we obtain for the neutrino mixing matrix 



R> = -0.52 -0.42 -0.74 (42) 




2 v 2 2 

with eigenvalues m 1 ~ 0.007-^-, m 2 ~ 0.12-^- and mi ~ 1.02-^|-, for M 3 the heaviest right- 
handed neutrino Majorana eigenvalue. The small contributions to the mixing angles coming 
from the charged lepton sector, Eq. (jHZJ) and Eq. ()44j) . are only significant for the CHOOZ angle 

sinfl 13 ~ j sin0^ - sintfg s[ 2 e iipL | (43) 
From here we see that the CHOOZ angle can be reduced or enhanced depending on the phase 

Regarding CP violation in the neutrino sector, we have three different observables, the MNS 
phase measurable in neutrino oscillations, the leptogenesis phase and the neutrinoless double 
beta decay phase. 

In first place let us determine the MNS phase. First we define an hermitian matrix = 
Sfc XikX'jt with x' the low energy effective neutrino mass matrix defined in Appendix|Bl Eq. (JHSj) . 
We calculate the MNS phase through the invariant Im (H" 2 H 23 H 31 ) in the basis in which the 
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charged leptons are diagonal 9 . Then we notice that we can factor out the phases P%Pj 1 from 
H?j and these phases cancel exactly in the invariant. Then the only remaining phase up to order 

e 2 - 75 is ifL = —^O-d — 2a u + (3 3 . In this way, we can calculate the invariant from x" = P~ 1 x'P~ 1 i 
which from Eq. (|88jl can be written at leading order in s 12 , as 

/ Xu ~ 2 X 2is[ 2 e i ^ X12 - (X22e i<PL - Xue' iipL )s[ 2 Xis ~ XvA^ \ 
x"-[ - X22 + 2 X i2S l 12 e- i ^ X23 + Xiss'ue-^ , (44) 

\ ••• ••• X33 / 

for s\ 2 ~ eX\/T, 2 . Then the H'X elements (as computed from x") can be written as follows 

H» = A l2 + s[ 2 e^{\ X ii\ 2 -\X22\ 2 + \Xn\ 2 -\x^\ 2 ) 
H'» = A 31 -s[ 2 e- l ^A 32 

H» = A 23 + s[ 2 e-^A 13 (45) 

where we have defined Ay = J2 a =i 3 XiaX*j a an d we must remember that in this case Ay are real 
up to order e 2 - 75 . Then 

H V 12 H^H» X « A 12 A 31 A 23 + s[ 2 [A 12 (A 2 31 - A 2 3 )e"^ + 

A 3 iA 23 (|xn| 2 - \X22\ 2 + \xrs\ 2 ~ |X3 3 | 2 )e^] + . . . (46) 

where there can be higher powers in s l l2 e ±lipL , although it is enough to keep the powers shown. 
Due to the hierarchies in the elements Eq. (JHEJ), we have that xu « X22, X23, X32, X33- Then 
we can approximate, 

MH^H^H^) ~ Im(Ai 2 A23A 3 i) + [^12X23 (X22X33(X22 - X33) + xlaixm + Xas)) + 

X13X23X22 (X2 3 (X33 - X22) + X22X33 (%22 + X33))] s[ 2 Sm (f L (47) 

This expression is completely general and includes the effect of non-negligible charged lepton 
mixings. In fact, we have A12A23A31 = 0(x 2 a Xbc) with a, b, c = 2, 3 while the square braket in 
Eq. (|4T|) is 0(xiaXbc)- Therefore, it follows that the charged lepton phase will always contribute 
at leading order if s l 12 > XiaJ Xab which in fact is usually satisfied in a model where charged lepton 
and down quark yukawas are related. With the exception of the CHOOZ angle, the neutrino 
mixing angles are not affected by charged lepton mixing angles. It is important to emphasize 
that there is no equivalent result when considering leptonic phases |2~?t I28j. As we see above, a 
relatively small charged lepton mixing is enough to affect at leading order the MNS phase. 

In our particular case, taking into account that A 12 A 2 3A 31 is real and replacing the elements 
of the effective neutrino mass matrix, Eq. (|86|h in the case 77 = r/, we obtain, 

MH^HM) « s 1 ^ sin *l (l - £ c 2 2 X 2 ) 4 ^ s 12 4(1 + s 2 2 ) X?X 2 (X 2 + X 2 ) (48) 

Provided, as is generally true, this imaginary part does not vanish the observable MNS phase is 
given by this term, = — 2«d — 2a u + (3 3 . 

9 Note that after some simple algebra, using the unitarity of the MNS mixing matrix, Im (H'^H^H^) — 
«-<)(<-<)«-<) Im (V2 2 V 2 * 3 V 33 V 3 * 2 ). 
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Table 2: Transformation of the superfields under the SU (3) family and R x Z-i x U(X) symmetries which restrict 
the form of the mass matrices. 



The leptogenesis CP asymmetry is of the form 

*i = — — r — - — y Im [(Y^Y^U ^ (49) 

where the neutrino Yukawa matrices are given in the basis of diagonal Majorana masses, Eq. ()85|). 
In this case, it is easy to see that the diagonal phase matrix cancels; in particular the phases 
coming from charged leptons do not contribute. Then the phase relevant for leptogenesis is given 
by if = Aa u + f3 3 . 

Finally, the neutrinoless double beta decay phase is simply the relative phase between the 
two dominant contributions of x'n i n t ne basis of diagonal charged lepton masses. This time 
from Eq. (jHHj) we obtain, 

(X')n = Pi (xn - 2f|W^) (50) 

and if we remember that Xij are rea l U P t° order e 2 ' 75 we obtain the result that the neutrinoless 
double beta decay phase is also v?l i.e. it coincides with the MNS phase. 



4.2 Model II 



As a second example we will use the model proposed in reference This model is based in a 
gauged SU(3) x 5*0(10) broken to a Pati-Salam subgroup. The additional symmetries needed to 
obtain the correct Yukawa textures and Majorana masses are given by a Z 2 x U(l) x R symmetry 
with the charges specified in Table El With these symmetries the Yukawa superpotential is, 

w Yuk ~ [0\0{ + 0^0{^ + (e M e 23tk e 3il 0i 3 + e jk % 3 , k e 3 ,tei 3 ) {0 23 W 3 ) + 

+ e^ k er 3k (9 23 e 3 ) 2 + ^% k (0 23 3 )(0 23 23 ) + (eyi + ei4 3 ) s] (51) 



Maj 



>VA' + l 23 J 23 (X0 23 )(X0 3 )(0 3 23 ) 3 + 



M 
1 



^^23^3, fc )(^ m #23^3, m )(A#23)(A# 3 )(#3#23)(#23#3 



(52) 



Now, if we compare the Yukawa superpotential, Eq. ()51|) . with the superpotential given in 
Eq (JHB*|) for Model I it is clear that the Yukawa superpotentials are indeed very similar and the 
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only difference is in the term that mixes O^O^. Clearly the magnitude is such that it does not 
affect the quark masses and mixing angles much, although we have seen it can be significant in 
determining the prediction of the phase $ 2 . On the other hand, the Majorana superpotentials 
Eq. PU|) and Eq. (J52j) are very different and this modifies the structure of neutrino masses and 
mixing angles. 

The vev of the A fields break lepton number and gives rise to the Majorana matrix, 



-iSx/2 \ / £ 6 f 3 q o \ / e ~ iSl/2 

e-^l 2 • e 6 e 2 ■ e"^/ 2 I M 3 . (53) 



e 6 e 3 










e 6 e 2 


i) 










a u + 


4/3 3 . 


Here 



as in the up sector due to SU(4:)ps and we have factorised two diagonal phase matrices that, as 
shown below, play no role in leptonic CP violation. 

Apply the see-saw formula again we obtain the low energy neutrino mass matrix. As in [Sj, 
using the analytic formulas in f2S] and [27], we obtain 



mi ~ e± (54) 

m2 „ k±|)!ji.4.8#- (55) 



2 



M 3 M 3 

Mi .U< 



m 3 « — — — ~ 13.5—— (56) 



tan 



tan Q\ 2 



®13 



1.2 



(g+j) 
(5-1) 



-cs 23 



0.6 



+ |)Kg ~ |)(-f ) + (9 + |)(-f + ce- 1 / 2 )] 
[G?-f) 2 + G? + f) 2 ] 3/2 



2e 



(57) 

(58) 
(59) 
0.05, 



where the numerical estimates correspond to g ~ 1, h ~ 0.3, c' ~ — c ~ 1, a ~ e ~ 
£ ~ 0.15. 

To obtain the CP violating phases in the neutrino sector we proceed as in the previous 
case. The effective neutrino matrix in the flavour basis (before diagonalising the charged lepton 
Yukawa matrix) is, 



X = Pu 



X l + e i(S2-h) XiX 2 



e 

X1X2 

e 



+ e 



+ e 



i(S2-Si) 



*(<5a— h) 
i(S2—Si) 



•Pu 



(60) 



/ 



with 5 X = Si + 25 u and P v = Diag(e i(<5 " +(52 - <5l/2) , e iSl/2 , e i{5l/2+M ). Here, we have neglected 
sub dominant terms suppresed at least by e 3 . Now we go to the basis of diagonal charged lepton 



17 




masses using Eq. ()87|) and we obtain, 

1 -ege*** \ / 1 
x ' = P- | £^e~^ L 1 • x • -£j£ eiVL 1 ] • P (61) 
e o 1 / \ 6 

where now (p L = 2(a u - a d ) + 4a u + 4/3 3 and P = ^^^(^-2^+8^/2)^^/2^ e i(s 1 /2-2 X +2p 3 )y 
In this basis we calculate the MNS phase through the invariant Im( H" 2 H 23 Hi 3 ). We can use the 
formula in Eq. (f4T)|) although, in this case the are no longer real and have a phase at order 
e with respect to the real part. In this case, = ^ a =i 3 XiaXja are 5 

A 12 = 2^ + X 1 (2 + X 1 )e- i(52 ^ l) + ... 
A 13 = 2^^ + X 1 (2 + X 1 )e i ^ ) + ... 

A 23 = ^^( X l + xl)+2^e^-^ + 2^e-^-^ + ... (62) 

e 2 e e 

However, we can easily see that Im(A 12 A23A 31 ) is 0(1/ e 3 ) while the leading contribution to the 
invariant is, 

« 4^^p(X? + X^)(X* - X^ + 2X?X^) sin^ L + . . . (63) 

Zj e £ 

and the real part of the invariant is given at leading order by AXfX^Xf + X|)/e 4 . Therefore 
in this case the MNS phase is given by ipL = 2a d — 2a u — 2j3 3 + 2(5' 2 and again charged lepton 
phases play a dominant role in the MNS phase 

Now, we must calculate the leptogenesis phase from Eq. (J^)) in the basis of diagonal Majorana 
masses, Eq. f!85|) . Again in this case, the phases coming from charged leptons cancel and therefore 
this phase is unrelated to the MNS phase. 

Finally the neutrinoless double beta decay phase is the relative phase between the two dom- 
inant contributions to x'n- We have 

(X')n = p i (xn " 2s^- Xl2 e-^ ~ PlXle^^ (l - 2 ± e ^-^ (64) 

Therefore, the neutrinoless double beta decay phase is simply, (S 2 — $1 — <Pl) = 2(a u — a d ). This 
is just the CKM phase in the quark sector. 

In summary, the phases are constrained in specific models although the constraints are model 
dependent. In Model I the MNS phase, measureable in neutrino oscillations, is the same as the 
phase relevant to double beta decay processes. However there is no relation with the quark CKM 
phase in this case. In Model II the double beta decay phase is just the CKM phase but now 
is not related to the MNS phase. In both cases the leptogenesis phase is not given in terms of 
the other phases. As we shall discuss in Section® all the phases are determined at the stage of 
spontaneous breaking. This leads to a discrete set of possible values for the phases. For example 
in Model I the CKM phase is quantised in units of 7r/2 showing that maximal CP violation in 
the quark sector is quite natural. All the other phases are similarly quantised but the unit of 
quantisation is too small to make phenomenologically interesting predictions. 
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5 Supersymmetric soft breaking and the flavour and CP 
problems 

So far we have seen that an SU(3) model of flavour with spontaneous CP violation can success- 
fully reproduce the quark and lepton masses and mixing angles together with the observed CP 
violation. We have also calculated the predictions for CP violation in the neutrino sector. In 
this Section we analyse the new effects of these phases in the supersymmetric sector. 

5.1 The SUSY CP problem 

Perhaps the most severe problem in supersymmetry phenomenology is the so-called "supersym- 
metric CP problem" . This refers to the SUSY contributions to the electron and neutron electric 
dipole moments (EDM) that in the presence of phases in the \x term or the trilinear couplings 
typically exceed the experimental bounds by two orders of magnitude. Therefore, first we must 
check these phases. 

The \x term in our model is generated through a Giudice-Masiero mechanism. This means 
that there is no bilinear in the Higgs fields in the original superpotential and this term is only 
present in the Kahler potential as XHiH 2 + h.c. Once supersymmetry is broken an effective 
/i term is generated in the superpotential of the form Xm^^HxHi- In our model, CP is an 
exact symmetry at Mpi an k, before the breaking of the flavour symmetry. Thus, at this scale, the 
term /i = Am 3 / 2 is real jS]. Now we have to worry about the corrections proportional to flavon 
vevs once we break the flavour symmetry. However, it is easy to see that these corrections will 
be extremely small and irrelevant for phenomenology. First we must take into account that in 
the absence of the /i term (or equivalently of the mixed Higgs term in the Kahler potential) 
our theory has a new U(l) symmetry distinguishing the two Higgs doublets. Therefore any 
correction to the \x term must contain the \i term itself. This implies that the first correction to 
the \x term must come from a two loop diagram [32] ■ Then only a neutral combination of flavon 
fields under both SU(3) and the other flavour symmetries (different from the trivial ones 9 a 9l) 
can contribute and we have seen in the previous Section that the size of this neutral combination 
is constrained from the structure of the Majorana mass matrix and can be at most e 2,75 . In 
summary, the phase of the \i term is well below 10~ 4 . 

Regarding the trilinear couplings, we can obtain the trilinear couplings after the flavour 
symmetry breaking in terms of the effective superpotential Eq. (J2J) and the effective Kahler 
potential Eq. (|5jl. Here, it is much more convenient to work in the basis of diagonal Kahler 
metric for the visible fields. In this basis, the trilinear couplings, {Yp)ijHpQiq c ^ are, 

(Tp)a = yfjF a d a (K + \og{KlKlK])) + F a d a Yf 3 (65) 

where is the Yukawa matrix in terms of the canonically normalised fields, Yf- are the original 
Yukawa couplings before cannonical normalisation and the Kahler potential is K = K(X, X') + 
K\{X, X^)\(j) l \ 2 with X hidden sector fields, visible sector fields and K\(X,X^) real. 

The first term in Eq. (|65|) gives rise to a factorisable contribution to the trilinear couplings 
Yff = (Ai + Aj)Yj EH with At = F a d a \ogK\ and Aj = F a d a (K + ]og(K$K%)) real. From [HI, 
and taking into account that these Ai are real and the only phases are in the Yukawa matrix, 
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we can see that the contribution from these factorisable nonuniversal trilinear terms to EDMs is 
safely below the experimental bounds. The remaining term in Eq. (|65|) involves the derivative 
of the Yukawa couplings in terms of fields with nonvanishing F-terms. If we make the derivative 
in terms of the flavon fields themselves it was shown in [HT] that the diagonal trilinear couplings 
in the SCKM basis are real at leading order in the flavon fields. In the same spirit, we can prove 
that this is also true for derivatives in terms of real hidden sector fields. If the Yukawa matrix is 
Y(9, X) = V L (9, X) ■ D(9, X) ■ V R (9, X) this contribution to the trilinear couplngs in the SCKM 
basis is, 

F\Vl ■ d a Y ■ V R ) U = F a {V ] L ■ o ViO«A + F a d a D, t + F a A(<9Vj ■ V R ) U (66) 

Now a given eigenvalue or an element of the mixing matrix can also be expanded in terms of 
the small flavon fields, Vy = J2 n (a n + b n X)9 n . If we expand Eq. (j^J in term of 9, the order in 
9 of the eigenvalue will set the leading term of the expansion as it appears in the three terms in 
Eq. (|66|). In fact, if A = (a m + b m X) 9 m for some m, d a Di will keep the same order in 9, oc 9 m . 
Moreover, as the only phases are in the flavon fields, d a D% has the same phase as A- Then the 
first or third terms in this equation, will give a leading order contribution to the expansion only 
if (V^ ■ d a V) does not depend on 9 and and in this case is completely real. So, all the leading 
order contributions keep the same phase as the eigenvalue itself and therefore is rephased away 
when we make the mass real. Only subdominant terms in the 9 expansion can get any phase, 
exactly as in |31| . This suppression is already enough to satisfy the EDM bounds. 

5.2 The SUSY flavour problem 

We have shown that the SUSY CP problem is elegantly solved in this theory. However, we must 
also check the status of the SUSY flavour problem (including offdiagonal CP violation). 

In this case we have two possible sources of new flavour structures associated with the 
breaking of the flavour symmetry, either D-term or F-term contributions. First, we must consider 
the possibility of large D-term contributions that can spoil the required degeneracy among 
squarks. In our models, the various D-terms associated to SU (3) fi breaking are small due to 
the symmetry of the flavon superpotential in the fields 9^,-theta^ and 6*23-^23 |5. 3 as shown in 
section El So we can neglect these D-term contributions to the sfermion masses. On other hand 
we also have contributions to soft masses due to some nonvanishing F-terms, as we have already 
discussed in Section 3, Eq. (JTTJ). 

Thus, we obtain the sfermion mass matrices from the Kahler potential Eq. (jHJ), as Mfj = 
m l/2 + F a F b (d log Kij) I (dcf) a d(f) b ) . The general structure of these matrices is given by Eq. (fTTj) 
although the expansion parameter is different for different SU(2) R quantum numbers e ~ 0.15 
for the right handed down quarks and charged leptons and e ~ 0.05 for right handed up quarks 
as well as all SU{2)l doublets (for a Pati-Salam group below Mpi an k). This means that the 
largest effects will be present in the right-handed down squark mass matrix and the right- 
handed charged slepton mass matrix, although differences in the left handed mass matrices 
can be also significant. At this point, we should emphasize that measurement of the complete 
flavour structure of the sfermion mass matrices will be the final check of the nature of the flavour 
symmetry. Ideally we would like to produce directly the different sfermions and simply measure 
their masses and mixing angles. However, before we are able to do that at LHC or at the Next 
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Linear Collider we can still look for indirect signs of SUSY flavour in FCNC or CP violation 
experiments [2]. 

Using Eqs. (j5J) and we can obtain the right-handed down squark mass matrix in the 
SCKM basis where we can use the phenomenological Mass Insertion bounds. However, this 
is not enough in the case of complex contributions, as we have still the freedom to rephase 
simultaneously left-handed and right-handed quarks without changing the eigenvalues. This 
rephasing does change the phase in different elements of the CKM mixing matrix and in fact 
in the presence of non universality new Yukawa phases become observable [33] ■ The MI bounds 
calculated in [20] are only valid in a basis where the phases in the CKM matrix are located 
in the (1,3) and (3, 1) elements (up to order e 4 ) 34j. Therefore we must go to this particular 
basis to compare our predictions with the bounds. More explicitly, the down Yukawa matrix in 

diag 



Eq. (0) is diagonalised, Y d = V[ ■ Y diag ■ where, 
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(67) 



and 
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diag R 



■v, 



where the diagonal phase matrices clearly do not affect the Yukawa matrix in Y = Vj~ ■ Y^ 
but combined with different rephasings in the up sector bring the CKM mixing matrix to the 

standard form with phases in (1, 3) and (3, 1) elements and tu = Arg ( 1 



Now we have (M% ) SCKM 
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(69) 



In this case, we obtain a large {5rr)i2 from the rotation of the "large" difference among (1, 1) 
and (2, 2) elements, that generates the order e 3 entry in the (1, 2) element. The modulo of this 
entry is (for e = 0.15 and Xi = 1.4) is 5 x 10~ 3 at Mqut- For the imaginary part we obtain 



Im(M?J 12 



,d\2 



e^- 



-Xisin^xmp ~ 1.15 x 10 



-3 



rrin 



(70) 



where we taken $1 = 90° and |af| 2 /S u ~ 1.5. We must take into account that off diagonal 
elements of squark mass matrices do not change much from Mqut to M\y- However the diagonal 
elements get a large gaugino contribution and then, at the electroweak scale we can roughly 
approximate the average squark mass by 5mg. So we finally estimate the value of the Mass 
Insertion at the electroweak scale as, Re(5^ fl )i 2 ~ 1 x 10~ 3 and Im(5^ R )i 2 ~ 2.3 x 10~ 4 . 
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This must be compared with the phenomenological MI bounds [211] • In this case, the relevant 
bound is, 



\ model / \ bound 

y/lMSfudtt) = 6 - 8 >< 10_4 < (yiM^)? 2 lJ = 3.2 x lO" 3 (71) 

In a similar way, we can obtain the left handed down squark mass matrix in the SCKM basis. 
We must start from this matrix in the flavour symmetry basis, that now has a smaller expansion 
parameter, e ~ 0.05. 




M 2 3l ~ 1 + e 2 e 2 e 1 ^ m 2 (72) 
\ e 2 e~ i ^ 1 + £ J 

and then with the rotation Eq. (JHEJ) we obtain in the SCKM basis, {M 2 f)j ) SCKM = ■ M? • V[, 
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\n d \ 2 \n d \ 2 \ 
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2->d 

\n d \ 2 

y -ee 2 e- t( - 2x -^ l -^ L X 1 e V 2 * - eeY, d 1 + e j 



m 2 (73) 



The presence of offdiagonal entries in these squark mass matrices can have sizeable effects 
in several low energy observables. We have already analysed the effects of RR mass insertion 
in 6k and we have seen they can reach a 20%. However, if we consider simultaneously the 
RR and LL mass insertions the contribution could be even larger. For this combination the 
phenomenological bound is, A/|Im(5^, i? ) 12 (5^ L )i2| < 2.2 x 10~ 4 , while in our case even with the 
smaller LL mass insertion we obtain at the electroweak scale, 

|Im(^ R ) 12 (^ L ) 12 | = Vim ((1 + % 0.23)10- 3 (1 + i 0.23)10~ 4 ) = 2 x 10" 4 (74) 

for $i = 90°. This can give a large part of the experimental result. However, we must take 
into account that this represents the maximum possible contributions and all the elements in 
the sfermion mass matrices are multiplied by unknown factors order 1. Nevertheless we can 
conclude that a sizeable contribution to €k can be expected. 

In the case of $13 MI that could contribute to B^B^ mixing and the J/ipKs CP asymme- 
try, we can see from Eq. (J69j) and Eq. (J73|) that these MI are exactly of the same order as 
the corresponding 6u MI. However, the values for the MI required to saturate these observ- 
ables are now much larger and the phenomenological bounds are, y/|Im(<5^, R )i 3 | < 0.3 and 
a/| ^(5^)13(5^) 12 1 < 9 x 10~ 3 Thus no sizeable effects are possible here. 

Finally we can consider b — > s transitions with a S23, contributing to B S B S mixing and some 
CP asymmetries as B — > QK$ jSE]- From Eqs. (j69|) and ()73)1 we obtain at the electroweak scale, 

(5 d RR ) 23 ~ (4.5 + 6.8^) x 10- 3 (5 d LL )23^ (1.5 + 0.5 z) x 10- 3 (75) 

where, once more, we assumed maximal phases. Although in this case there are no MI bounds 
as in [201 from the literature on SUSY contributions to this decay EH], it is clear that 
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Field 
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Table 3: Transformation of the superfields associated with the spontaneous symmetry breaking. 



an isolated LL or RR mass insertions needs to be at least 0(0.2) to have a sizeable effect. It is 
not clear in the case of simultaneous LL and RR MI, but from [SHI ^ seems that a MI of order 
0.05 is required. Therefore, apparently there is no sizeable effect here. 

Lepton flavour violating decays, as \i — > and r — > /17, can also receive sizeable contribu- 
tions. In this case, the slepton mass matrices are exactly identical to Eqs. (jp^j) and (fT3~j). with the 
only replacement of — > S e (in this case the presence of phases is not important). However, 
the main advantage of leptonic processes is that the MI are not largely reduced from M GUT to 
My/-, as the diagonal elements of slepton mass matrices are less afected by gaugino masses in the 
evolution between these two scales, in fact at Mw we can approximate the average left handed 
slepton mass by 2mg. In this case, the bounds on the leptonic MI are more difficult to obtain 
as they depend on other parameters (//, tan/3, Mi/ 2 , ml). However, there some bounds for fixed 
values of tan/5 |3*D] . The most sensitive precess is fi — > ej where for tan/3 = 10 and average 
slepton mass of 300GeV we get (<5£x)i2 < 3 x 10~ 4 while for the RR mass insertion the bound 
is much worse due to a possible cancelation among diagrams. Now we have for the LL MI from 
Eq. (J73J), (S1 L )i2 < 2.5 x 10~ 4 and therefore a /1 — > decay close to the experimental bound is 
indeed possible. 

Finally, we have to comment on possible nonuniversality effects in the trilinear couplings. 
Trilinear couplings are obtained from Eq. (J65|) and in this case, it is easy to see (101 that 
large effects can be expected again in the kaon and leptonic sectors due to the big experimental 
sensitivity in these experiments. However these effects are usually sufficiently small in the B 
system due to the associated suppression with quark masses and the smaller sensitivity in B 
experiments. Again, as shown in [SI], these effects can be close to the experimental bound in 
the decay fi — > e'j if the flavon F-terms are sizeable. 

6 Vacuum Alignment 

It is easy to construct a superpotential capable of generating the phenomenologically desirable 
pattern of symmetry breaking vevs. Since in [S] a suitable superpotential for Model II has been 
constructed, we concentrate here on constructing a suitable superpotential for Model I. We will 
use it to demonstrate how maximal CP violation may readily be obtained. 

There is considerable freedom in constructing the breaking sector because it is sensitive to 
additional family and SM gauge singlet fields. Such fields abound in compactified string theories 
so we make no apologies for adding several such fields to the theory. An example of a set (we 
have not searched for the minimal set) giving rise to acceptable spontaneous breakdown is given 
in Table El The resulting superpotential has the form 
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W = P{9 3 9 3 + T) + U ({6 23 e 23 ) + S 2 ) + V ((9 3 9 3 ) 4 + SS) + 
+ Y(6 3 9 2 ) + VS 2 {9 3 9 3 ) 4 + Z {(9 23 9 3 )(9 23 9 2 ) + S 2 ) 
+ fxH 1 H 2 [1 + (9 23 9 3 )(9 3 9 3 ) 3 + T 5 + ...] (76) 

The overall symmetry breaking here is triggered by the fields 9 3 , 9 3 which obtain a vevs along a 
flat direction through radiatively breaking due to the Yukawa coupling P9 3 9 3 in Eq. ()76j) . Then 
it is transmitted to the rest of the flavon sector through the other terms in the superpotential. 
Note that the fields 9 3 , 9 3 couple symmetrically so that, if their soft masses are degenerate at the 
Planck scale, they will remain degenerate. As we show below, this is important in ensuring the 
close equality of their vevs which is necessary in avoiding large family dependent D-terms and 
the associated flavour changing processes [5]. Notice however, that a possible non- degeneracy of 
the soft masses would only change the relative size of the vevs in factors 0(1) but never their 
order in e and our predicted Yukawa textures would not be affected. 

From this flavon superpotential, we have the following scalar potential, 

V = \F V \ 2 + \F V \ 2 + \F V \ 2 + \F P \ 2 + |F y | 2 + \F Z \ 2 

+ m 3/2 /i E X E 2 [1 + {9 23 9 3 ){9 3 9 3 f + T 5 ] (77) 

where the last line corresponds to the soft terms associated with Eq. ()76)1 . Note that we have 
included in Eqs. f!76|) and (J77J) a /i term which only appears after SUSY breaking through a 
Giudice-Masiero mechanism. Notice that although this contribution to V is extremely small 
it is the first contribution lifting the flat direction left by the F-terms. Clearly, any neutral 
combination of vevs can multiply the superpotential but the terms in Eq. (|7B|) are enough to 
specify all the relevant phases in this case. 

After the fields #3, 9 3 get radiatively induced vevs, minimisation of Fy generates vevs for S 
and 5* which are equal, (5*) = (S), if they have equal soft breaking masses. Then minimisation 
of F v forces the alignment (9 3 9^)/M 2 = [9f9f)/M^ and minimisation of Fjj generates a vev for 
($23^23) • Using this together with D-flatness, we obtain, 

\(W _ KM 
M u M u 

1(023)1 _ I (fta) I 

M d M d 

Minimisation of Fz forces 9 2 to acquire a vev while minimisation of Fy makes it orthogonal to 
9 3 , by definition the 2 direction. Thus we see that minimisation of Fz forces 9 23 to get vevs in 
the 2 and 3 directions. Finally, inclusion in the scalar potential of the SUSY breaking soft mass 
terms for 9 23 (due to the ST/ (3) family symmetry they are equal in the 2 and 3 directions) shows 
that 9 23 will get equal vevs in the 2 and 3 directions - the required vacuum alignment. 

A major problem in models with a continuous gauged family symmetry is that the associated 
D-terms could split the degeneracy among different sfermions. Consider, for instance, the fields 
9 23 and 9 23 . Making Fjj zero fixes the vev of the product (9 23 9 23 ). Then, if the soft masses of these 
fields are equal we can also cancel the D-terms choosing equal vevs for the barred and unbarred 
fields. If the soft masses are not equal the size of the D-terms are roughly determined by the 
difference of the soft masses. Therefore, to minimise the value of the D-terms, we must assure 
that the soft masses of the relevant fields are reasonably close at the scale of symmetry breaking. 
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(78) 
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These soft masses must be approximately equal at M P i ank and this degeneracy must not be 
spoiled by radiative corrections between M P i ank and the symmetry breaking scale. Regarding 
the required level of degeneracy, -D 23 is precisely the most dangerous D-term, 

lAd 2 = 9 2 \0l^l; - ^3,2^2*3,3 + MSI', (79) 

as this would give an unsuppressed off-diagonal contribution to the squark masses. However, 
as seen in the previous section, the phenomenological bounds on this off-diagonal entry are not 
too strong. Moreover due to large gaugino contributions in the running from Mpi ank to My/ to 
diagonal squark masses, the relative size of the off-diagonal entries is reduced. Therefore, we 
only have to require a mild degeneracy O(0. 1) to the 823 and 823 soft masses 10 . This level of 
degeneracy is quite natural in string theory where groups of fields are expected to be degenerate. 
Another possibility that can also readily occur in string theory is that the soft masses of 823 
and 823 are < 0(0.1) m~, with nig the soft mass associated with the squark fields. Even if we 
assume that these masses are exactly equal at Mpi ank the degeneracy will be broken by radiative 
corrections if these fields have different interactions. However, as we have noted above, the fact 
that that the superpotential in Eq. (J76)) is approximately symmetric in 83 and 83 or 823 and 823 
implies their soft masses will not receive significantly different radiative corrections from Mpi ank 
to the symmetry breaking scale. As a result [5] the SU(3)fi D-terms are small and consistent 
with the bounds from FCNCs. 

Finally, following the discussion of [5], the lepton number violating vevs of A and A are 
automatically aligned along the 3 direction. 

We turn now to the really new feature, namely the spontaneous determination of the phases. 
They are also fixed from the minimization of the scalar potential. We first remind the reader 
that, due to the underlying CP symmetry, the suppressed couplings in Eq. (|77jl are real. Thus 
the minimisation of F terms involving the cancellation of two contributions requires that the 
relative phases of the two terms must be ±1. Then from the minimisation of the F-terms in 
Eq. (|77jl we obtain 



Further from the soft terms in the scalar potential, 

/3 3 + Aa u = 171x71 
3(a u - a d - x) + (a« - oc d ) = m 2 rr 

5 8 T = m 3 vr (81) 



with rii, rrii integers. 



8 T « 
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®s + 8 S = 

A(a u -a d -x) = 
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a u + m ir 
P' 2 + n 2 7i 
4 a u + n 3 it 
n 4 Ti 

a u + (3 3 + 5 2 + n 5 n 



(80) 



10 This D-term would also generate an offdiagonal entry in the (1,2) element of the squark mass matrix after 
diagonalisation of the corresponding Yukawa matrix. However, the phenomenological constraints are satisfied 
with the same level of degeneracy. 
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From this system of equations we obtain, 

7T 3 

T = 777,3- a u - a d = (m 2 - 7 n 4 )7r, (82) 
4 

The CP violating phase in the CKM matrix is $1 = 2(a u — a^) = —\m^K. We see that it is 
quantised in units of n/2 and so it is quite natural to obtain maximal CP violation in the quark 
sector in good agreement with observation. 

One may proceed in this manner to determine the other phases. Unfortunately, although all 
the other phases are also quantised, the unit of quantisation is too small to make phenomeno- 
logically interesting predictions. 

7 Conclusions 

In this paper we have investigated spontaneous CP violation in models based on an SO(10) x 
SU(3) fi symmetry. We have shown that such models provide a remarkably consistent description 
of the known masses and mixing angles of quarks and leptons. In the quark sector, the presence 
of CP violating phases is necessary, not only to reproduce CP violation processes, but also to 
reproduce the observed masses and mixings. We have also shown that the spontaneous breaking 
of CP in the flavour sector naturally solves the supersymmetric CP problem and the SUSY 
flavour problem, although flavour changing processes must occur at a level close to current 
experimental bounds. 

We have presented two possible variations of the model which simultaneously reproduce the 
observed neutrino masses and mixings. In these models there are phenomenologically inter- 
esting relations between the CP violating phases in the quark and lepton sectors and testing 
these relations will be a sensitive discriminator of models. The magnitude of CP violation is 
determined by the potential driving the spontaneous family symmetry breaking. The structure 
of this potential is such that it can quite naturally lead to maximal violation in the quark sector, 
consistent with present observations. An important complication, that we expect to be quite 
generally true, is that charged lepton phases, even if there is relatively small charged lepton 
mixing, cannot be ignored in the determination of the MNS and neutrinoless double beta decay 
phases. On the other hand the leptogenesis phase is independent of the charged lepton phases. 

The main prediction of the models based on a non-Abelian family group is the full flavour 
and CP structure of the sfermion mass matrices. Accurate measurement of sfermion masses 
and mixing angles will provide a significant test of such models. Even before we are able to 
access the sfermion sector at the LHC or at the Next Linear Collider there are several FCNC or 
CP violation experiments receive which sizeable SUSY contributions in such models and which 
already provide an experimental probe of the underlying family symmetry. 
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A Experimental Constraints 

We have performed a fit to the CKM elements (A, X, p and ff) with the experimental information 
available in Summer 2003, as appear in Table 0] In Table 1 of you can find other parameters, 
entering the formulas in the fit, which are well measured and we take as fixed in the fit. The 
relevant formulas for this fit can also be found in reference Figures ^ and 121 show the results 



Fitted Parameters 2003 


Parameter 


Value 


Gaussian-Flat errors. 


Refer en. 


A 


0.834 ±0.036 




* 


X 


0.2240 ±0.0036 






\v ub \ 


37.35 x 10" 4 


(±1.69 ±5.0) x 10~ 4 




\Vcb\ 


41.1 x 10~ 3 


(1.63 ±0.8) x 10- 3 




B K 


0.86 


0.06 ±0.14 




m* 


(1.3 ±0.1) GeV 




* 


m* t 


(167 ±5) GeV 




* 


rfx 


1.38 ±0.53 
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0.47 ±0.04 




E2 


Am Bd 


(0.502 ±0.006) ps- 1 




HH] 


fs d \/ B Bd 


0.228 GeV 


(0.030 ±0.020) GeV 


t m 




1.21 


0.04 ±0.05 




Am Bs > 


14.4 ps" 1 at 95% C.L. 




[43 



Table 4: Fitted Parameters 2003. The parameters marked with * have been computed as in [S] with the 
new data from |2T]. In fg d \jBs d (t) we have conservatively doubled the theoretical (flat) error from 021 to 
accomodate possible large results from unquenched lattice QCD calculations U2| 



of the fits in the (p, ff) plane. In all the plots, we have used the constraints from | V^/V^l, Am^ 
(lower limit), Am^ and |e#|. The constraint from sin 2(3 has not been included in the fit. 

In Figure HI due to the possible supersymmetric contributions to ex as presented in Section 
O we have enlarged the experimental error of a 20 % of its mean value, e#. Then, the 
SM contribution lies within the range ex ± 0.20^. We have considered that ex has a uniform 
probability within this range. However, we do not find a significant variation of p and fj, with 
respect to the original SM fit. Hence the phase $2 does not change considerably in this case. 

In Figure |2l we have assumed a fixed sign of the SUSY contributions, such that the mean 
of the experimental ex receives a SUSY contribution of 20% of its value, thus decreasing the 
SM contribution to 80% of ex- In this case, the fitted value of fj decreases (fj = 0.306^q!o52) 
with respect to the value (fj = 0.328^o'g3g) obtained in the original SM fit, allowing for the lower 
values of $ 2 (= 0.599lH)- 
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P 

Figure 1 : Fits of the parameters p, fj to the 2003 values of the constraints afecting the CKM matrix in the SM 
and in the MSSM case. The dotted lines indicate the regions of la for the different constriants, except in the case 
of sin 2(3 where we have ploted the 2a region as a reference. In \ejc\ the red(light) lines represent this constraint 
in the SM without new physics contributions. The grey(dark) lines represent the constraint in the MSSM case 
where we enlarge the error a 20% as explained in the text. The coloured (different grey tone) regions are the CL 
at 99%, 95% and 68% in the SM fit, while the wide continuous lines are the same CL in the MSSM case. 
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Figure 2: In this figure we assume that 20% of the measured value of efc is due to SUSY contributions and we 
fit only the remaining 80% of e^. As we can see in this case we can reach lower values of fj, required for lower 
values of $2- 



B Neutrino mass matrix, Model I 

With the seesaw mechanism, we obtain the low energy effective neutrino mass matrix \ v as > 

r = Y u Mft 1 (Y U ) T = Y»L M (M R )^ g L T M (Yn T , (83) 

where Lm is the complex orthogonal matrix diagonalising Mr and up to 0(e 5 ' 5 ), 

ic 12 e- i( - 5 '- 4au) s 12 e- i{5 '~ 4au) 
,. w - i -is 12 e- i4a " c 12 e- i4a " ^e 2 - 11 >e~ i4a - \ (84) 
isles' 2 - 75 -curjxe 2 - 75 1 
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where s 12 ~ and c 12 ~ 1 - ^e 1 / 2 . Now (M fl )^ g = Diag.(^, ^nr, 1) and 

( e i(5«-4a„) 

-zs 12 £ 3 (X! - X 2 'q 1 £ 2 - 75 e iv ) cx 2 e 3 (*i - X 2 r ]l s 2 - 75 e i *) e 3 (X 2 e ilp + X^e 2 - 75 ) 

ie 3 (ci 2 X 1 + s 12 rirne 2 - 5 e^) e 3 (si 2 Xi - ci 2 ^i£ 2 - 5 e^) r/e 2 - 75 e^ | , (85) 

ic 12 s 3 X 2 - is l2 (r] - ri^e 2 - 75 s l2 e 3 X 2 + c l2 (r] - ij^e 2 - 75 1 



3- 



where we have defined p = (r] 2 — r/f), X\ — (1 + X), X 2 — (1 — X) and the phase 9? = 4a u + $ 
Notice that the structure of the rotated Yukawa matrix has allowed us to extract a diagonal 
matrix of phases such that only ip appears in the remaining matrix. Now we can obtain the ele- 
ments of the effective neutrino mass matrix vn v LL = v 2 P-x u ■ P with P = Diag(e^ <5u_4a "' ) , e~ l/3s , 1) 
and, 

X v ^ (86) 

( -M2X1 f,s 12 c l2 X 2 {l + s 2 2 ) ^s 12 X, (c 12 X 2 (l + s 2 2 ) - { -^f^s\ 2 ) \ 

-i?X 2 c\ 2 -£ (x.X.c 2 , - ^c 12 (l + s^X,) 

V 1 + 2^c 12 (l + s 2 2 )X 2 - %c 2 2 X 2 - ^s 2 2 ) 

However this is the neutrino mass matrix before diagonalisation of charged lepton Yukawas. 
To obtain the lepton mixings and especially the CP violating phases relevant for phenomenology, 
it is more convenient to go to the basis of real and diagonal charged lepton masses. To reach 
this basis we must rotate and rephase this effective neutrino mass matrix: x' = L * • P • x ■ P ■ L , 
with, 



e 



IS* \ / ! g x 



L l = P' ■ U = I 1 • -efi 1 I + 0{e 2 ) (87) 

e i(ax-A) / V E 

which is enough for our purposes. So, we have, 

1 -£j±e iVL \ / 1 e Je _iVi 

x ' = P . I ege-*** 1 • x • -^fte i¥,i " 1 | • P (88) 

e o 1 / \ 6 

where <p L = -2a d - 2a u + (3 3 and P = Diag(e i (- 2a <- 2a ">, e~ il33 , e'Hfc+fc)) 
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